GMRT discovery of PSR J1544+4937, an eclipsing black-widow 
pulsar identified with a Fermi LAT source 
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ing black-widow millisecond pulsar, PSR J1544+4937, identified with the un- 
cataloged 7-ray source Fermi J1544. 2+4941. This 2.16 ms pulsar is in a 2.9 
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ABSTRACT 

Using the Giant Metrewave Radio Telescope (GMRT) we performed deep ob- 
servations to search for radio pulsations in the directions of unidentified Fermi 
Large Area Telescope (LAT) 7-ray sources. We report the discovery of an eclips- 
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hours compact circular orbit with a very low-mass companion (M c > O.O17M ). 
At 322 MHz this pulsar is found to be eclipsing for 13% of its orbit, whereas at 
607 MHz the pulsar is detected throughout the low-frequency eclipse phase. Vari- 
ations in the eclipse ingress phase are observed, indicating a clumpy and variable 
eclipsing medium. Moreover, additional short-duration absorption events are 
observed around the eclipse boundaries. Using the radio timing ephemeris we 
were able to detect 7-ray pulsations from this pulsar, confirming it as the source 
powering the 7-ray emission. 

Subject headings: pulsars: general; binaries: eclipsing, pulsars: individual (PSR 
J1544+4937) 



1. Introduction 



The Large Area Telescope (LAT; Atwood et al. 2009) aboard the Fermi Gamma-ray 
Space Telescope has discovered a large number of 7-ray point sources, of which many are 



unidentified or even unassociated with any known potential counterpart (Ackermann et al. 



2012). The LAT can localize most of these sources well enough that they can be covered in a 



single pointing with the large primary beams of radio telescopes at low-frequencies, allowing 
them to be searched efficiently. Targeted radio searches of unassociated LAT point sources 
by the Fermi Pulsar Search Consortium (PSC) have resulted in the discovery of 43 radio 
millisecond pulsars (MSPs; Ray et al.|2012 ). MSPs are thought to evolve from normal pulsars 
in binary systems via transfer of angular momentum from companions. Thus, the majority 
of MSPs are naturally expected to be in binaries (~ 83% being the binary fraction for MSPs 
in the Galactic fieldj. Binary systems where the pulsar wind evaporates the companion 
are one way to form isolated MSPs. Such systems where the interaction is ongoing are 
called black-widow (BW) pulsars. Many exhibit long eclipses (~ 10% of the orbital period, 
apparently larger than the companion's Roche lobe) that are believed to be caused by the 
material blown from the very low mass companion (M c <C O.1M ) by the pulsar wind. There 
were two such eclipsing BW systems in the Galactic field known before the launch of Fermi 



PSR B1957+20 ( jFruchter et al.||1988[ ) and PSR J2051-0827 ( |Stappers et al.||1996[ ). The 



BW pulsars are found to have higher values of spin-down energy-loss rate (E ~ 10 34 erg s x ) 
compared to other MSPs, making these systems good candidates for pulsed 7-ray emission 



(Roberts 2011). Among 43 new MSPs found in Fermi-directed searches there are at least 



1 http: //astro. phys .wvu.edu/GalacticMSPs/ 



10 BWs (Ray et al. 2012). This Letter describes the discovery and follow-up study of an 



eclipsing BW MSP, J1544+4937, with the GMRT. 



2. Observations and search analysis 

As a part of the PSC search effort, we observed mid- and high-Galactic- latitude unas- 
sociated Fermi point sources with the GMRT at 607 MHz. The GMRT Software Back-end 
(GSB; |Roy et ah 2010) produces simultaneous incoherent and coherent beam filter-bank 
outputs of 512x0.0651 MHz sampled every 61.44 /j,s. The wider incoherent beam of the 
GMRT (beam-width 40' at 607 MHz) can easily cover error circles associated with the Fermi 
sources. In addition a coherent beam that is 3x more sensitive and narrower (1.5' at 607 
MHz) can be useful if the pulsar happens to be near the pointing center. 

One of the targets was Fermi J1544. 2+4941, a 7-ray source from an unpublished 
internal source list created by the LAT Collaboration using 18 months of data in preparation 



for the 2FGL catalog (Nolan et al. 2012). The source location (J2000) from that analysis 
(and used for our telescope pointing) was R.A. = 236°. 074, Decl. = 49°. 695, with a 95% 
confidence error circle of radius 9.5'. This source is very weak, with a test statistic (TS) of 



26.2 (Mattox et al. 1996) in the 18 month analysis, and did not make the significance cut to 



be included in the 2FGL catalog itself. 

We processed the data on an IUCAA HPC cluster with Fourier-based acceleration search 



methods using search techniques implemented in PRESTO (Ransom et al. 2002). We inves- 
tigated trial dispersion measures (DMs) ranging from pc cm -3 up to 350 pc cm -3 . A linear 
drift of up to 200 Fourier-frequency bins for the highest summed harmonic was allowed. The 
powerline, 50 Hz, and its subsequent harmonics were excised. Using parameters of 32 MHz 
bandwidth, 10% duty-cycle, incoherent array gain of 2.3 K/Jy, for 30 minutes of observing, 
we estimate the search sensitivity as (92K +T s k y )/(335K) mJy for a 5cr detection at 607 
MHz. Considering |6| > 5°, where T s k y ~ 10-45 K, our search sensitivity is 0.3-0.4 mJy. 

In a 30-minute pointing on 2011 February 1, towards Fermi J1544. 2+4941 we dis- 
covered a binary MSP of period 2.16 ms with significant acceleration of 2.25 m s~ 2 at a DM 
of 23.2 pc cm" 3 . 



3. Follow-up timing 

We localized J1544+4937 with an accuracy of 5" (positions listed in Table ffl) using 
continuum imaging followed by multi-pixel beamforming (Roy et aL]|2012), which allowed us 



to have sensitive follow-up studies using the GMRT coherent array. We estimate a flux of 
5.4 mJy at 322 MHz, and a spectral index of —2.3. We started the regular timing campaign 
for Jl 544+4937 in April 2011 at 322 MHz. With the derived position from the multi-pixel 



search and an a priori binary model predicted by Bhattacharyya & Nityananda (2008) 



we obtained phase-connected time-of-arrivals (TOAs) from TEMPOrl using the JPL DE405 
solar system ephemeris (Standish||2004). The binary timing model used is ELL1 (Lange et al. 



2001), since J1544+4937 is in a very low eccentricity system. This MSP is in a very compact 



binary with an orbital period of 2.9 hours. We derive a minimum companion mass (for 90° 
orbital inclination) of 0.017 M Q using the Keplerian mass function, assuming a pulsar mass 
of 1.4 M . J1544+4937 is eclipsed for about 13% of the orbit at 322 MHz (SecQ. The 
best-fit timing model (MJD 55680.927-56332.90) is obtained excluding the TOAs around 
the eclipse phase. We achieved a post-fit rms timing residual of 6.9 /xs from 652 days of 
timing (Fig. QJ). We estimate a precise DM equal to 23.2258(1) pc cm~ 3 by a timing fit using 
322 and 607 MHz TOAs from non-eclipsing binary phases (excluding TOAs between phases 
0.05—0.35). The resulting ephemeris, position from timing fit and derived parameters are 
listed in Table [Q 
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Fig. 1. — Post-fit timing residuals of J1544+4937 considering non-eclipsing binary phases. 
The best-fit timing solution is given in Table [1} 



Table 1. Parameters of J1544+4937 



Parameter Value a 

Interferometric position 13 

Right ascension (J2000) 15 E 4^)4?llxmJ?3 

Declination (J2000) +49°37'57'.'45±4?7 

Offset from survey beam center 4.3' 

Parameters from radio timing 

Right ascension (J2000) 15 h 44 m 04?48722 (2) 

Declination (J2000) +49°37'55'.'2545 (2) 

Position Epoch (MJD) 51544.0 

Pulsar period P (ms) 2.15928839043289 (5) 

Pulsar frequency / (Hz) 463.11553585462 (1) 

Frequency derivative / (Hz s _1 ) -6.29 (l)xlO -16 

Period Epoch (MJD) 56007.0 

Dispersion measure DM (cm -3 pc) 23.2258 (1) 

Binary model ELL1 

Orbital period P b (days) 0.1207729895 (1) 

Projected semi-major axis x (lt-s) 0.0328680 (4) 

Epoch of ascending node passage T A SC (MJD) 56124.7701121 (2) 

Span of timing data (MJD) 652 

Number of TOAs 280 

Post-fit residual rms (fis) 6.9 

Reduced chi-square 2.7 

Derived parameters 

Mass function / (M ) 0.0000026132 

Min companion Mass m c (Mq) 0.017 

DM distance (kpc) 1.2 

Flux density at 322 MHz (mjy) 5.4 

Flux density at 607 MHz (mjy) 1.2 

Spectral index —2.3 

Surface magnetic field B 3 (10 8 G) 0.805 (1) 

Spin down luminosity E (10 34 erg s" 1 ) 1.150 (8) 

Characteristic age r (Gyr) 11.65 (3) 

7-ray parameters 

Photon flux (> 0.1 GeV, cm" 2 s" 1 ) 1.6(8) xlO -9 

Energy flux (> 0.1 GeV, erg cm" 2 s" 1 ) 2.1(6) xl0~ 12 

Luminosity, L 7 //n (> 0.1 GeV, erg cm -2 s" 1 ) 3.6 xlO 32 

Efficiency, jj 7 // n (> 0.1 GeV) 0.03 

a Errors in the last digit are in parentheses. 
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Fig. 2. — Phase-aligned radio and LAT 7-ray light curves of J1544+4937. The horizontal 
dashed lines are an estimate of the background level from sources other than the pulsar. 

4. 7-ray pulsations 

The radio timing position of J1544+4937 is 4.3' from the LAT localization of Fermi 
J1544. 2+4941. This is well within the radius of the 95% confidence error circle, suggesting 
an association between the pulsar and the LAT source, which we sought to confirm with 
detection of pulsations. In addition, the 7-ray detectability metric E 1 / 2 /^ 2 for J1544+4937 
is 7.4 x 10 16 erg 1 / 2 kpc -2 s -1 / 2 , which is comparable to other 7-ray detected MSPs (Fig. 12 



of Abdoet al. 2010). 



Because of the low significance of the LAT source, a spectral analysis of the region 



http : //www . atnf . csiro . au/research/pulsar/tempo 



is required to optimize the sensitivity of the H-test for pulsed significance using photon 
probability weighting (Kerr 2011[ ). We performed a binned likelihood spectral analysis of 
the region using the Fermi Science Toohjj We selected LAT data with reconstructed energies 
in the range 100 MeV to 20 GeV collected between 2008 August 4 and 2013 February 6. 
We used the P7SOURCE_V6 instrument response functions and excluded events at a zenith 
angle > 100° and times when the LAT was in the SAA or the rocking angle was greater 
than 52°. We modelecfj a region of radius 15° using a source list from the 2FGL catalog 
with an additional source added at the location of the pulsar. We modeled the pulsar with 
an exponentially cutoff power-law spectrum. Because it is so faint, we fitted only for the 
power-law index and normalization, keeping the cutoff energy fixed at 2.5 GeV, a value 



typical of other MSPs (Abdo et al. 2010). The spectral model included both the isotropic 



and Galactic diffuse contributions with normalizations free in the fit. We detected the source 
with a TS of 38.2, corresponding to ~ 6er. The spectral index and cutoff energy are very 
poorly constrained for this weak source, so to estimate the uncertainty on the integrated flux 
we repeated the spectral fit with the cutoff frozen at 0.7, 2.5, and 5.0 GeV, and with the 
index frozen at 1.3 (the average value for LAT-detected MSPs) and the cutoff free. We find 
that the uncertainty due to the spectral shape is about 3 x 10~ 10 cm~ 2 s^ 1 for the photon 
flux and 2 x 10~ 13 erg cm -2 s _1 for the energy flux. These are smaller than the statistical 
errors for this faint source and are added in quadrature to compute the errors on the fluxes 
reported in Table [I] Using the fitted energy flux (G), and the DM distance (d), we compute 
the 7-ray luminosity L 1 = Airf^Gd 2 and efficiency r] = L^/E, where the beaming factor /q 



is assumed to be 1 (Waters et al. 2009). 



Using this model, we selected events within 2° of the source and computed the proba- 
bility that they originated from the source with the Fermi Science Tool gtsrcprob. Using 
these probabilities as weights and pulse phases computed using the radio timing model, we 
computed a weighted H-test of 37.1, corresponding to a detection significance of 5.1 a and 
confirming that J1544+4937 is indeed a 7-ray pulsar and identifying the LAT source as this 
MSP. The LAT weighted lightcurves in two energy bands, phase aligned to the radio profiles, 
are shown in Fig. [2] For this faint source, it is difficult to determine the peak multiplicity, 
but the 7-ray emission seems to be mainly between phases 0.4 and 0.6, relative to the single 
radio peak. Fitting a single gaussian peak to the light curve gives a peak at phase 0.54 ±0.02 
with a width of 0.25 ± 0.05. 



http: //f ermi . gsf c .nasa.gov/ssc/data/ainalysis/software/ 



http : //f ermi . gsf c . nasa . gov/ssc/data/access/lat/BackgroundModels . html 



5. X-ray and optical observations 

We analyzed 5 Swift/XKT observations (totaling 9.3 ks) of this field from 2013 March 
19-24. No X-ray counterpart is detected, providing a flux upper limit of ~ 1.7xl0~ 14 ergs 
cm" 2 s" 1 for 0.3-10 keV. 

We have checked for an optical counterpart for J1544+4937 using archival SDSS data. 



This field was imaged (Adelman-McCarthy et al. 2008) on MJD 52046.294, at binary phase 



0.57. No optical counterpart was detected. The strongest upper limits were g' > 22.5, 
r' > 22.0 and i' > 21.5 (90% CL). We also obtained a 300s Ha image with the MiniMo 
camera on the 3.5m WIYN telescope on MJD 55975.521 (Brownsberger fc Romani|2013). No 



stellar counterpart or extended emission was seen with an Ha flux upper limit of 2.5 x 10 16 
erg cm -2 s -1 (90% CL) for a compact < 5" nebula. 

We compare these limits to fluxes measured for the archetypal, and best studied, BW 
B1957+20 ( |Reynolds et al.|2007|[v"an Kerkwijk et al.|2011[ ) which has a DM distance ~2.5kpc 



(with an estimated extinction Ay = 1.2), P& = 9.17 hours and E = 7.7 x 10 34 ergs -1 . At 
the orbital phase corresponding to J1544+4937's upper limits (0.32 past eclipse) the studies 
above give extinction corrected magnitudes g' ~20.2, r' ~19.7 and i' ~ 19.4. Assuming a 
Roche lobe filling factor similar to that of the B1957+20 (since we do not know the true 
size of the companion), the ratio of reprocessed heating flux is just the ratio of spin-down 
luminosities, i.e. J1544+4937 is 6.4 x lower. Whereas, if we bring B1957+20 to 1.2 kpc, it 
would be 4.3 x brighter. The extinction correction for J1544+4937 is very small (Ay = 0.05) 
in this direction. Corrected magnitudes of B1957+20, scaled to J1544+4937, would be 
0.4mag fainter. Observed limits on J1544+4937 are Ag' = +1.8mag, Ar' = +1.9mag and 
Ai' = +1.7mag, fainter than the scaled B1957+20 values. Thus the optical data imply 
~ 2.3 x larger distance of J1544+4937. Alternatively, the companion radius may be smaller 
by a similar factor or the companion temperature may be lower than the heating comparison 
above. 

Our Ha limit for J1544+4937 corresponds to < 0.15 of the highest surface brightness 
~ 5" patch of the B 1957+20 bow shock. The average flux ratio is even more constraining, 
with an upper limit of ~ 0.03 of the total bow-shock flux. However, given the small fraction 
of MSPs that show optical bow shocks (likely due to the small filling factor of the neutral 
interstellar medium), the non-detection at such a large distance from the Galactic plane (b~ 
50°) is not particularly surprising. 



6. Eclipse characteristics 

Fig. [3] presents the timing residuals (and electron column densities) around the eclipse 
phase at 322 and 607 MHz, with 90s time resolution. The effect of the eclipses is generally 
seen from 0.18 to 0.31 orbital phase (eclipse-zone hereafter) at 322 MHz. The eclipses are 
centered at binary phase 0.24 with a duration of around 22 minutes. We estimate the radius 



of the companion's Roche lobe, Rl (Eggleton 1983) 



49ao 2 / 3 
fli = 0.6^» + ]■(!+,■/») ~ °- 1Me (1) 

where q = m c /m p is the mass ratio of the companion and the pulsar, and a is the separation 
of the companion from the pulsar (a ~ 1.2i? Q for J1544+4937, indicating an extremely 
compact binary). The opaque portion of the companion's orbit is 0.98 R , much larger 
than Rl of 0.13 R , so the volume occupied by the eclipsing material is well outside the 
companion's Roche lobe, and thus is not gravitationally bound to the companion. This 
confirms that this binary is a BW, where the pulsar is ablating its companion, creating a 
significant amount of intrabinary material that obscures the pulsar's emission. 

Our sample consisted of six eclipses at 322 MHz where the pulsar emission was fully 
obscured by companion and its wind and two 607 MHz observing sessions covering the 
full orbit (Fig. g. At 607 MHz we detect the MSP throughout the 322 MHz eclipse-zone. 
However, we observe a flux fading at 607 MHz near the superior conjunction (orbital phase 
0.24). Significant delays in pulse arrival times are observed at 607 MHz during the eclipse- 
zone at 322 MHz. Maximum delay in pulse arrival time at 607 MHz near the eclipse superior 
conjunction is around 300 /is, which corresponds to an increase in DM of 0.027 pc cm -3 and 
an added electron density, N e of 8xl0 16 cm -2 . 



Thompson et al. (1994) (T94 hereafter) elucidate a collection of eclipse mechanisms. 
According to them, eclipsing due to refraction of the radio beam demands an order of mag- 
nitude higher group delay (~ few tens of milliseconds) than we observe for J1544+4937 (250 
/is at 322 MHz egress). Using N e ~8xl0 16 cm -2 observed during superior conjunction at 
607 MHz and absorption length L about twice the size of the eclipse-zone, according to 
equation 11 of T94 we find that free-free absorption is possible (absorption optical depth 
Tff > 1) if the plasma temperature T< 4 x fd ' K, where f c i ' is the clumping factor. This 
demands either a very low temperature or a very high value of the clumping factor, both 
of which are not physically achievable. Eclipsing by pulse smearing (due to the increase of 
N e along the line of sight) can be ruled out, as the excess electron column density inferred 
from 607 MHz predicts 374 /xs smearing of pulses at 322 MHz near superior conjunction 
(considering incoherent dedispersion), which is less than one fifth of pulse period. Since 
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J 1544+4937 has a narrow main-pulse, and no significant profile evolution is apparent at the 
eclipse boundary, pulse broadening due to scattering can reduce the detectability but cannot 
explain the eclipse. In addition, since J1544+4937 is relatively weak, nearby and has a shal- 
lower spectrum than B1957+20, the expected induced Compton scattering optical depth is 
much less than one (equation 26 of T94). Another eclipse mechanism considered by T94 is 
cyclotron-synchrotron absorption of the radio waves by non-relativistic/relativistic electrons, 
which requires a magnetic field in the vicinity of the companion. We calculate a magnetic 
field B ~ 11 G and corresponding cyclotron absorption frequency ~ 31 MHz (equations 35, 
37 of T94, assuming a moment-of-inertia of 10 45 g cm 2 ). Thus 322 and 607 MHz will corre- 
spond to 10th and 20th harmonics of the cyclotron resonance. For a fixed temperature, the 
optical depth for cyclotron absorption drops with harmonics, which may explain the lack of 



absorption seen at 607 MHz. A modified model is proposed by Khechinashvili et al. (2000) 
based on kinematic treatment of cyclotron damping, assuming white-dwarf companions with 
reasonably strong surface magnetic fields. Further observations over a wider radio spectrum 
may help to investigate the frequency-dependent degree of damping predicted by this model. 

We observe a temporal variation of the ingress phase at 322 MHz. For two eclipses, the 
ingress phase is shifted to 0.16 (from 0.18 for the other four eclipses), whereas there are no 
apparent shifts in the egress phase. This corresponds to an increase of the opaque portion 
by 0.15 R Q and N e >1.2xl0 16 cm -2 at the ingress boundary. Such asymmetric increase of 
eclipse duration is unique to this BW and may indicate that our line of sight is probing a 
wind zone where there is systematic outflow of eclipse material. 

We also observe strong phase modulations and additional short-duration absorptions at 
ingress and egress, in time-series data at higher resolution. The durations of these features 
are in general around 10-20 s, and hence they are not seen in Fig. [3} However, in one 
of the observing epochs these modulations lasted longer — phase modulation of duration 
100 s, followed by a short- duration absorption of ~ 180 s, then regular emission resumes 
for 500 s, after which the eclipse starts (Fig. [5]). Such intriguing modulations, indicating 
fluctuation of plasma density around eclipsing zone, are not reported for other eclipsing 
BWs. Fragmented blobs of plasma randomly oriented around eclipsing zone, obscuring 
radiation from the pulsar, can explain these short-duration absorptions. 



7. Discussion 

We report the GMRT discovery of an eclipsing BW pulsar, at the position of an unasso- 
ciated LAT source, Fermi J1544. 2+4941. This is the first Galactic field MSP discovered 
at the GMRT. The detection of pulsed 7-rays from this pulsar demonstrates it as the source 
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powering Fermi J1544. 2+4941. Due to the limited significance of the source in 7-rays 
additional data are needed before conclusions on peak multiplicity and system geometry 
can be drawn. The implied efficiency (~ 3%) of converting spin-down energy into 7-rays is 
typical of LAT-detected MSPs. This is the first discovery of a radio MSP in a LAT source 



fainter than the 2FGL catalog limit (Ray et al. 2012|. Since the radio pulsar is relatively 



bright, this provides strong justification to continue these searches as new LAT sources are 
revealed in analyses of longer datasets. The radio flux is uncorrelated with the 7-ray flux 



(Ackermann et al. 2012), so even faint new LAT sources can harbor bright radio MSPs. 



By assuming a comparable companion Roche lobe filling factor and orbital inclination, 
we find that the irradiated flux/unit area should be 0.73 x that of B1957+20, leading to an 
effective temperature of the heated side of T max ~ 7700 K. The flux limits from section [5] 
imply a distance > 2.7 kpc. This would imply a larger 7-ray efficiency of > 15%, which is 
still typical of LAT-detected MSPs. Deeper optical observations with multicolor light curves 
of the companion can be used to test these estimates. 

Eclipsing BW pulsars have the potential of providing information on the evolutionary 
connection between the low-mass X-ray binaries and isolated MSPs. With long term mon- 
itoring of this pulsar we aim to estimate Pb and its higher derivatives, which can provide 



an estimate of the life span of the system. Bates et al. (2011) noted that for BW systems 
the measured value of E/a? is an order of magnitude higher than for other MSP binaries, 
indicative of greater energy flux needed to ablate the companion. For J1544+4937 we cal- 
culate E/a 2 ~ 1.5 xlO 33 erg lt-s _2 s _1 , which is similar to other BW systems in the Galactic 
field. Dual frequency observations presented in this paper suggest that cyclotron absorption 
by the plasma formed via interaction of the pulsar wind with ablated material can obscure 
the pulsed emission. However, exploring the radio spectrum on either side to probe the re- 
duced/increased opaqueness of the stellar wind during the eclipse phase may provide better 
insight into the plausible eclipse mechanism. 
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Fig. 3. — Variation of timing residuals and electron column density (TOAs of 90 s time 
resolution) around eclipse phase at 322 MHz (top) and 607 MHz (bottom). 
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Fig. 4. — Example of additional short-duration absorption beyond the eclipse phase at 322 
MHz on 2012 June 3. Left: variation of timing residuals and electron column density at 42 
s time resolution. Right: phaseogram showing short-duration absorption. 



